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Introduction 


In a previous paper [1] we have discussed variational methods for many-body 
systems with strong two-particle interaction with the use of wave functions built 
up from two-particle wave functions. General variational equations were derived, 
and simplified methods with parameter variation were discussed. One-dimensional 
systems were studied in some details, and for a system of bosons with Gaussian 
interaction explicit expressions were calculated with a simple choice of trial wave 
function. 

As a preparation for the study of more realistic problems, we now make a more 
detailed analysis of a one-dimensional system of fermions with Gaussian interac- 
tion and with the corresponding simple trial wave function, in order to examine 
how the method works with strong short range interaction and to use the one- 
dimensional case as a starting point for the extension of the method to the 
three-dimensional case. Explicit calculations are made for a system of four three- 
dimensional fermions, and we close our work with a discussion of the application 
of the same method to the nuclear case. 

When working with trial wave functions for many-body systems that are not 
built up of Slater-determinants of single particle wave functions, one usually has 
troubles to avoid complicated integrals. Even for the simplest trial wave func- 
tions our integrals will contain lots of terms, and in order to carry through the 
integrations without using approximate methods whose correctness is usually hard 
to estimate, it seems to be important to find integration methods that make use 
of all the symmetries in the expressions. Some details of our computations are 
given in an appendix. 


System of one-dimensional fermions 


We consider a system of » one-dimensional fermions with Gaussian interaction 
potential and a common oscillator potential. The Hamiltonian of the system can 
then be be written 


hr S a.m iis tate > ex {—B (a,—a,)* (1) 
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If we use the “Overlapping pair method” [1] and choose as two-particle wave 
function the function describing two particles with oscillator interaction and a 
common oscillator potential : 


f (a, %_) = (%, — xg) exp {— kg (m+ tq)” — ko (4 - t)"}; 
then the trial wave function for the whole system is 
Hy Aye ety) 
TH (ey) ere Baden) 


xexp {—(n—1)(k@ + ko) aw — (ke —ke) ) 2, MM) 


where kg and ko are the variational parameters that will be determined from 
the equation 


6<H) =0. 


The calculation of the explicit expression of (H> = E (kg, ko) was made in (I) with- 
out much labour for the boson case, but in the fermion case the integrals will 
be complicated, due to the fact that the n-dimensional orthogonal transformation 
that transforms the exponentials to diagonal form must be applied to the poly- 
nomial A%, of the degree n(n—1), and this gives a great number of terms even 
for small n. For details of the integration, see the appendix. The expression for 
E as a function of the new dimensionless constants and parameters x, o, y and 
v, defined by 


2(n—I(k@tko)=x——", — 2(k@—ka)=7——, =dvho] 
nf dears and x—v=0 | 


Saye ine lyse 
jee Le o ies o ) 
2 oty “\oty/]]’ 
where P,, is a polynomial. In the appendix the coefficients are given for n= 2, 3 
and 4. 


This expression for the total energy of the system will now be minimized with 
respect to the parameters x and o. 


oH 
AEE gives nx—(n—1)o=1 or x+(n—1)r=1. 
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Fig. 1. The dependence of structure-determining parameter ¢ on the strength of the interaction 
v, is given for linear systems with two and four particles. If the figure is extended to higher 
v-values, the curves intersect again for v/y*15. The energies of the system H,, and HE p obtained 
by the variational method and by a first order perturbation calculation, are also given as a func- 
tion of v for the four-particle system. The line of short dashes indicates the smallest strength of 
the interaction v,,;,, for which bound two-particle states exists. A comparison of the two energy 
calculations shows that the first order perturbation can be good for an interaction strength up 
to4v but in the vicinity of v,,;,, it seems to be very bad. The interaction range is chosen 
rather short, y = 20. 


min’ 


The interpretation of this is that the center of gravity motion corresponds to an 
oscillator motion with the frequency m. We have then 


2] 2. =a. { 
E(c)=tho a (c+ )+2e IV. - Pa( zs 4 
2 Oo 2 ot+y o+y 


= ~ 


@EL/éc=0 gives the relation between v and o 
tk iy a 1/o? -1 en) 
n @ [V/ o P,( oO | 
da|! ory ot+y | 
Without interaction between the particles (v=0) we have o=1, and the wave 
function describes, as required, a system of independent Fermi-oscillators with 

the frequency w. 

The dependence of o on v is given in Fig. 1. o appears here as the single 


structure-describing parameter, and its deviation from 1 determines the deviation 
from independent particle state. 
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Disregarding the common oscillator potential, we get the expression for the 
total energy 


B= ho [Sat bn(n eee - Bal = )|. 
2 2 ot+y ory 


As function of % it has a minimum on the boundary ko =0, that is 
nu—(n—1)o=0, 


and from this follows that the wave function is independent of the center of 
gravity coordinate. The minimizing of this expression with respect to the parameter 
o gives the relation between v and o 


(n+ 1)/n 
a Ces ' 6 
relicerstitecd| 


and the final expression for the energy 
; =] é 
H=ho [boty Ma y V/ er | |. 
= 2 oty ot+y 


The least value of » for which solutions of o exist, can be interpreted as the 
weakest interaction for which bound states exist. 

Ordinary perturbation theory, with the two-particle interaction term as pertur- 
bation term, gives to the first order 


ial 1 
E—E,= th ee P. as) : 
owl lias 
This energy expression is compared with the corresponding expression obtained 
from the variational treatment in Fig. 1. : 


Three-dimensional system of Fermi-particles 


In the three-dimensional case it seems to be hard to find any two-particle 
trial wave function that gives a simple parameter wave function for the whole — 
system with the required property of giving the correct wave function for a sys- 
tem of free particles for certain values of the parameters, if we use the ‘“‘over- 
lapping pair model” for the building up of the system wave function. With the 
“disjoint pair model” (see part I) it is certainly possible to construct two-particle 
wave functions that give such a system wave function. This wave function, how- 
ever, properly antisymmetrised, consists of a sum of terms each lacking symmetry 
between the particles, and is therefore hard to handle. For the following discussion 
of three-dimensional systems we therefore choose to give up our starting-point that 
the wave function for the system will be constructed from two-particle wave 
functions, and instead choose the system wave function in direct analogy with 
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the linear case but retain the requirement on the wave function to describe an 
independent particle state for certain parameter values. 
We note that our system wave function in the linear case is of the form 


W= A, (x, one sted aE CP aves. 9 Uns Pi: Po: pets 


where A, is the antisymmetrical basic invariant [2] in the variables 2,, ..., 2. 
This is the antisymmetrical polynomial of the lowest degree in these variables. 
S, is a symmetrical function that contains the parameters p,, p., .... This is of 
course a general way of writing an antisymmetrical function, even if we restrict 
all functions to be regular. (For A, has only the zeros required by the antisym- 
metry.) For the three-dimensional case we now choose the corresponding form 
for our trial system wave function 


Be AL, com (Fal ion (lis Shins Pra Dar sae) > 


where A, is now antisymmetrical in the three-dimensional vectors 7; and a poly- 
nomial of the lowest possible degree in the vector-components, and S,, is a sym- 
metrical function. If we now restrict the functions to be regular, this is not a 
general way of writing a regular antisymmetrical function, but we have restricted 
our choice of trial wave functions to a class of functions with some prescribed 
ZeYOs. 


Three-dimensional system of four spinless fermions 


For four particles A,(r) is uniquely determined and rotational invariant. If we 
now choose as a simple example a system with central Gaussian interaction be- 
tween the particles and a common isotropic oscillator potential, we must assume 
the ground state to be an S-state and must then choose S, to be rotational in- 
variant, too. 

As a trial wave function we now use 


W=NA,(r)S, 
and in analogy with the linear case we choose 


S,=exp { = DEY, 


where 
n 

= 2 — 9 \e] = 2 : : 
bW)=X thatnttkgin—nl= > [eXatty > mms @) 

x and o are given by (2). 

The Hamiltonian of the system is 
— SA ot 3 oh du exp (Blin) (4) 
sn ecpie ign ixj a Wide 


| The calculation of the expectation value of the total pon) is now very easy 


_ (see appendix), and we get the following expression for H in the same para- 
Bracters as in the linear case: 
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ok 
? 


Fig. 2. The dependence of the parameter 6 and the energy H on v is given for the three-dimen- 
min* 


sional four-particle system with two different ranges of the interaction, y=5 and 20. When the 
See text to Fig. 1. 


interaction is of short range, y>5, there is a discontinuous change in 06 when v passes a value 
slightly higher than v 


it esi fs 3 
E (#,0)=tho [3(2x+0)—120(—2_) +4(2+ 
if 


a 4x- ale 
The minimizing of this expression with respect to x gives the relation 4x—30=1 
which means that the center of gravity of the system moves as an oscillator 
with the frequency w. We have then 


1 / \ 5/2 
rel) 
om ory 


and the second minimizing equation, é#/éo=0, gives the relation between v and o 


a oa 
4 oyVo/(o+y) 
Neglecting the common potential, we get 


1 (o+ y)8 


4 oy Va/(o+y) 


and this function has a minimum for. o/y 


sce BB 

=§, 
as the weakest potential for which bound states exist. 
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which may again be interpreted 
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The dependence of EH and o on v is given in Fig. 2. 


Particles with spin 


The straight forward extension of the method to particles with spin is to in- 
troduce the ordinary spin-coordinate o with the possible values +1 and —1, and 
to define the antisymmetrical basic invariant as the polynomial in the variables 
Xi, Yi, %, and o;,i=1,...,, that is antisymmetrical in index exchange and is of 
the lowest possible degree in the Cartesian coordinates. 

For a system with spin independent Hamiltonian we then make the correspond- 
ing assumption for the wave function 


Ramee Cassese) 2g APO) eT cial s, Kgs 3) x Dar x02) 


where S, is independent of the spin coordinate. For certain values of n, A, (r, c) 
is not uniquely determined, and then the ground state level is degenerated. A 
weak spin-dependent interaction will usually split this degeneracy, but we can 
assume that the main influence of the spin-dependent term on the ground state 
wave function is on the choice of the right A,, and that S, will be the same 
function for all these degenerated states, and that the parameter values can be de- 
termined from the Hamiltonian without the spin-dependent term. For strong spin- 
dependent forces it is certainly not possible to get a useful trial wave function 
in this manner, but a different A, must be chosen, and S, may perhaps depend 
on the spin coordinate. 


Suggestion for the nuclear problem 


To follow the same line in the construction of a trial wave function for the 
atomic nucleus we construct in the same way the basic invariant A, (7, o, tT) and 
try to find a suitable S, that depends on a number of parameters for the wave 
function 

V(r, 0, 7) =A; (7%) G; BT) Ss (1) Py) Dos ---)> 


It is now tempting to assume that S, is a function that varies in a uniform 
manner with the particle number and contains the main properties of nuclear 
matter (including also the size of the nuclei), and that A, that varies irregularly 
with n, describes the main part of the individual properties in the nuclear table. 
However, even if the nuclear forces are charge independent, there is certainly a 
strong spin dependence, and the coulomb force will have a strong influence on 
the wave function; so it is only for light nuclei that one can use the basic in- 
variants A, with any hope of success. If we retain our requirement of S, to be 
regular and of the wave function to describe an independent particle state for certain 
parameter values, we are forced to use other antisymmetrical polynomials than the 


basic invariants for heavier nuclei. Possibly one can construct these polynomials 


: 
: 
: 
7 


of higher degree than the lowest with the guidance of shell model calculations. 
If one uses such a nuclear wave function with a uniformly varying S,, it would 
be of interest to examine if the main results of the shell model calculations remain 


with an S-function which departs from one describing an independent particle state. 
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Comparison with the Jastrow method 


Jastrow [3] has proposed a variational approach to the quantum mechanical 
many-body problem, with a trial wave function for a fermion system of the form 


P= ITF (15) y gly 


where ’, is a Slater-determinant of single particle orbitals, and f(7;;) is a corre- 
lation function depending on the distance betweeen the particles 7 and 7. f(7;;) 
contains also variational parameters. Several authors [4,5] have used this form 
of trial wave functions for the study of the so-called nuclear matter. They use 
for the calculation of their integrals the so-called cluster development method, 
i.e. an expansion of the two-particle density function in terms containing inte- 
grals over an increasing number of f-functions. 

We can write our trial wave function for the linear system in the following 
form 

n n 
Y =exp (o—1) 2 wt+y > 2 xexp (> 2? | A; 


1<i<j<n 


and this wave function corresponds to a Jastrow wave function with 


fley)=exp (0-1) Sabt+e > mal, 
1 


1l<i<j<n 


and with the variational parameters o and ». 


APPENDIX 


Here are given some details from the integrations of the expressions for the 
expectation values of the energies. 


The linear case 


The integral for the energy-expectation-value is 


Baw? | An () exp {—}b(x)} H A, (x) exp {—}6(z)} dz, ...d%, 


where win | 4 exp {—b(z)} da, ... dap, 


NOE a em 


A, (x) = (Rb (a, — «;) a 


| be fe xe xe... ah 


and the function in the exponent 


b(z)=x> 2i+y> a,2, (ais  xJength). 
i i<j ho 
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The Hamiltonian H of the system is given by (1). The interaction term in the 
Hamiltonian gives an additional term, f > (#;—.2;)’, in the exponent, and we have 
i<j 


to solve integrals consisting of polynomials multiplied by exponential functions 
of the argument 


b’ (w, B)=x> Dore + 2 oi a, + B(x, — 2). 
t <j 


This quadratic form is now “‘diagonalised’”’ by the solution of the corresponding 
secular equation. The roots of this equation are simply obtained by elementary 
determinantal transformations. One gets the actin expression for the quadratic 
form in the new variables 


B(x, B) = (v4 (w= 1) BE + (x9 + 2) EF S GH — 


The direction of the new coordinate axis in the n-dimensional space can be rep- 


resented by the vectors é,, ...,@,, where é,=(1, 1, ..., 1) and the other vectors are 
orthogonal. This is sufficient to “diagonalise” b (x), but for b’ (x, 6) one must choose 
one of the vectors to have the direction (1,—1,0,...,0). A, can be written as 


the product of the n(n—1)/2 planes 
Ni; ‘T= 0 


with the normal vectors 7,;;= (0, ...,0,1,0,...,0,1,0, ...,0), where 1 and —1 stand 
on the coordinate places i and j. All the planes intersect in the line (1, 1, ..., 1). 
We have then 
A, = II (vi 


i<j 


We now get the expression for A, in the new coordinates, if we express the nor- 
mal vectors in their components in the new coordinate system with the direction 
WECUOLS, Cs x ure, Cn 

When @, and @, are fixed, one can choose the other direction vectors in dif- 
ferent manners. A general and simple way to construct them is 


Cre LaeietsL) Cc Sak aio '35 0,75. 6.30) 
alas, OF 0) Ce ei (Mise PLS ARO ity) 0) 
=(1,1,-2, ore cers) a.s.0. 
With these new axis directions we get the expression for A, in the new coor- 
dinates 
A,=[lay, 
i<j 
f=1 
where ay= —(i-1)&+ > Satis. 
ksit+ 


For certain n-values it is possible to choose the axis directions in a more sym- 
metric way. For n=4 we have used for our computations the axis directions 


é,=(1, By; 1), é =(1,—1, 0, 0), 
é, = (0, 0, Ua é,=(1, Det. — 1). 
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If n=2”, where m is an integer, it is possible to construct the vectors in a cor- 


responding manner. 
The integrations now give the energy expressions 


TY -¥ Ea (~—») oa +f, (m) +2 fs on) 


2 
and 
ms 1 2h | 
(Vet. = 4 Ee 1)» @—py) 
where 
fy (n) = nl [ 4232 exp (—D(E)}dE, ... dE, 
and re 


1 i] 
nat | 43, (6a) xn v0 e6 ats 


It is not necessary to carry out the integrals for f, and f,, when we observe that 
they are independent on x and o, and that with »y=1 and x=mw/h the wave 
function describes noninteracting Fermi oscillators, where we have <7’) = <Vpor.> = 
=ln?haw. We then get 


Rae orf ; 
NP not? 4 Oe in— lly (x=) 


The integral for the expectation value of the interaction energy gives 


Seoul ses Aes 
<V. » ) Py ’ 
(Vint. = Yo 9 a—v+B ao 


where P,, is a polynomial. A general formula for the coefficients of P, must be 
very complicated, and we have only calculated its coefficients for m up to four. 


If P, (x)= > ata’, 


we have for the coefficients a” 


DS 

\n 2 3 4 

v ox 

0 1 1.125 1.3032 
il — 0:75 — 1.9310 
2 0.625 2.4976 
3 — 2.4355 
4 1.5657 
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The three-dimensional case with four fermions 


The integral for the expectation value of the energy is here 


co 


1 
E= NE | A,(r) exp {—b(r)} H A, (r) exp {—b(r)} Br, d3r,d? 7,03 14, 
ly yy & 
eee eet) meee 
where A, Be 8 
a(t) 1 2% Ys 2 
Dig Ya Ry 


and 6(r) is given by (3) with n=4. H is given by (4). A,(r) can now be written 
A,(r)=c¢ > Ps, y,2(—1)? (fi, -&) (q+) (Ag * 2) 
=¢ > P1,2,3 (—1)? (hy -&) (fi J) (ig « 2). 


The sums are taken over all permutations P, and p is the number of inversions 
in the permutation. The bar here means four-dimensional vectors, and 7,, 75, fig 
and the axis (1,1, 1,1) are orthogonal. (The vectors 7 ought here not to be mutu- . 
ally orthogonal but linearly independent.) These directions are, properly chosen, 
also principal axis of the quadratic forms in the exponent, and we choose them 
as our new coordinate axis. We then get the simple expression for A, (r) 


ae = 2 Psy (= 1)? 25 Yg%q = > Pos. (— 1)? 22 Y3 24, 


where the numbers now mean vector components in our new four-dimensional 
coordinate system. 

Because of the simple invariance properties of the Hamiltonian for this coor- 
dinate transformation, the integrals are very simple, and we get without much 
labour 


2m 3 
mo" ( 15 3 
<V rot.) = 4 (==, eer 


The present investigation was carried out at the Institute of Theoretical Physics in Stock- 
holm, and I wish to thank Professor O. Klein, Director of the Institute, for his kind interest 
in this work. 


Institute of Theoretical Physics, Stockholm, November 1958. 
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